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Recent numerical simulations address a conjecture by Shapiro that when two neutron stars collide 
head-on from rest at infinity, sufficient thermal pressure may be generated to support the hot 
remnant in quasi-static equilibrium against collapse prior to neutrino cooling. The conjecture is 
meant to apply even when the total remnant mass exceeds the maximum mass of a cold neutron 
star. One set of simulations seems to corroborate the conjecture, while another, involving higher 
mass progenitors each very close to the maximum mass, does not. In both cases the total mass 
of the remnant exceeds the maximum mass. We point out numerical subtleties in performing such 
simulations when the progenitors are near the maximum mass; they can explain why the simulations 
might have difficulty assessing the conjecture in such high-mass cases. 



Shapiro jlj recently speculated that following the head- 
on collision of two neutron stars from rest at infinity, the 
two recoil shocks which propagate back into each star fol- 
lowing contact might generate sufficient thermal pressure 
to hold up the remnant against gravitational collapse, at 
least until it cools via neutrino emission. This quasi- 
cquilibrium state will last many dynamical timescales 
(~ msec ) because the neutrinos, which eventually carry 
off the thermal energy, leak out slowly, on the neutrino 
diffusion timescale (~ 10 sec). The argument was inde- 
pendent of the total mass of the progenitors and applies 
even if the mass of the remnant greatly exceeds the max- 
imum mass of a cold neutron star. A simple, analytic 
analysis using relativistic polytropes served to motivate 
the conjecture. The analysis assumed strict conserva- 
tion of rest-mass and total mass-energy in the collision, 
as well as relaxation, following merger, to the same poly- 
tropic density profile in the remnant as in the progenitors. 
Given these assumptions it was possible to show that a 
stable equilibrium solution always exists for the remnant, 
independent of mass. In evaluating these assumptions, 
Shapiro argued that while they greatly oversimplify the 
collision in order to permit an analytic proof, they do 
provide a reasonable approximate description of the ide- 
alized head-on scenario under consideration. He did show 
that the loss of energy due to gravitational radiation, 
though small, would rule out an equilibrium solution for 
configurations arbitrarily close to the maximum mass, if 
the collision were not accompanied by loss of rest-mass. 
More significantly, he cautioned that a dynamical colli- 
sion need not relax to a stable, equilibrium solution after 
all, even though one exists, but could instead overshoot 
the equilibrium state and collapse to black hole. But he 
put forward the possibility of delayed collapse as a plausi- 
ble outcome of head-on collisions from rest at infinity and 
one which might be realized in numerical simulations. 

Numerical simulations to test this conjecture have been 
reported recently in Ref The simulations treat col- 
liding r = 2 polytropes in a 3+1 numerical scheme that 
solves Einstein's field equations of general relativity cou- 
pled to the equations of relativistic hydrodynamics. The 



simulations employed Y = 2 polytropic models that sat- 
isfy the TOV equilibrium equations for isolated, spherical 
stars. The stars are placed at finite separation, no more 
than 3i? apart, where R is the stellar radius. These con- 
figurations are then boosted towards each other with the 
Newtonian freefall velocity ||. No attempt is made to 
correct for the distortions in the initial matter density or 
pressure profiles that arise at such close separation due 
to tidal interactions. The simulations described in Ref 
Q deal with collisions of 1.4M Q stars 0. We shall as- 
sume below that this quoted mass represents the total 
mass-energy of each of the progenitors, although no dis- 
tinction is made in Ref between rest-mass Mq and to- 
tal mass-energy M, nor is there any discussion about how 
the mass in the colliding system is actually measured and 
monitored numerically. The subsquent evolution of this 
model produces a black hole a few dynamical timescales 
following contact and merger. This example, it is claimed 
in Ref Q , is in contradiction to the conjecture of Shapiro. 

Two aspects of the numerical simulations are crucial to 
note. One is that the maximum stable mass of the poly- 
tropic equation of state is 1.46M Q in the units of Ref M. 
The other, alluded to but not reported in Ref O, is that 
a second set of simulations using the same equation of 
state for the head-on collision of 0.8M Q configurations 
did not result in collapse following merger Jj|. In this 
case the simulations show that the merged remnant is 
supported in stable equilibrium by the thermal pressure 
generated by collision-induced shock heating, in appar- 
ent accord with the Shapiro conjecture. Interestingly, in 
both cases the total mass of the merged remnant exceeds 
the maximum mass of a cold star. 

Does the simulation with 1.4Mq stars show unambigu- 
ously that a delayed collapse does not occur for stars of 
this mass? Insight into what could be happening may 
be gleaned from the original discussion in Ref [Q. Sec- 
tion VI of that paper, "Discussion and Caveats " already 
anticipated numerical subtleties that would be encoun- 
tered in simulating a collision of stars that are initially 
at rest at infinity. To obtain reliable results, the require- 
ments for computational accuracy become very stringent 
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for high-mass stars near the maximum mass. The rea- 
son is clear from the analytic argument for the existence 
of a stable, equilibrium solution for the remnant. For 
such a solution to exist, it is necessary that the bind- 
ing energy of the final remnant reside along the stable 
branch of the TOV equilibrium curve. The ratio M/Mq 
measures the specific binding energy of a star accord- 
ing to Ebmd/Mo = 1 — M/Mq. This ratio monotoni- 
cally decreases with increasing central density p c along 
the stable branch of the TOV curve, whereas the mass 
and rest-mass monotonically increase along this branch. 
The turning point on the curve marks the onset of radial 
instability; there M and Mq asssume their maximum val- 
ues and the binding energy ratio assumes its minimum 
value, (M/Mo) m i n . In the idealized scenario in which the 
stars begin from rest at infinity, strict conservation of 
mass and rest-mass imply that for the hot remnant 
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As a result, the binding energy of the hot remnant is 
identical to the binding energy of each of the progenitor 
stars, and this guarantees the existence of an equilibrium 
solution (cf. Eqs. 8 - 10 of Ref §). 

Now consider an alternative scenario in which Eq. (jl]) 
is replaced by 
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As discussed in Ref [0, the fractional change in the to- 
tal mass-energy / might be the result of a true physical 
departure from the conditions or assumptions of the ide- 
alized collision scenario. The quantity / may account, 
for example, for the energy loss due to gravitational ra- 
diation or for the release of the progenitor stars from rest 
at finite separation rather from infinity |6) . Alternatively, 
/ might represent the degree of numerical error arising 
in a numerical simulation. Such error might be present 
due to imprecise initial data, including, for example, ini- 
tial data which do not exactly correspond to the true 
solution at finite separation for two stars which actually 
begin from rest at infinity . Numerical error might also 
result from finite-difference integration error in the evo- 
lution due to finite spatial grid and time-step resolution, 
imprecise outer boundary conditions or outer boundary 
points that don't extend sufficiently far into the radiation 
zone, global instabilities, etc Q. To illustrate, suppose 
we assume that for whatever reason, numerical error re- 
sults in a spurious fractional decrease in M in the course 
of the simulation. Wc then estimate the the maximum 
tolerable error /toi > which still admits a stable, equi- 
librium solution for the remnant |}). According to Eq. 
(0) , relaxation to a stable equilibrium state requires 
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Eqn. (0) demonstrates the difficutly of simulating 
accurately a head-on collision from rest at infinity as 
the progenitor masses approach the maximum mass. As 
M -> M max , we have (M/M ) -> (M/Af )mm and, hence 
/toi -► 0. 

To appreciate the severity of the problem, consider 
the actual simulations performed for Y = 2, for which 
(M/M ) min = 0.910 (see Fig 3 in Ref @). This mini- 
mum occurs at 1.46M Q , the maximum mass. We then 
find (M/M ) = 0.915 for the M = 1.4M star and 
{M/M Q ) = 0.959 for the M = O.8M star. According 
to Eq. (||) these values imply 

/ to i(0.8M Q ) = 0.05 and / tol (1.4M ) = 0.005. (4) 

Eqn. (|j) suggests that determining whether stable equi- 
librium state forms following the merger of two 1.4Mq 
stars requires that the simulations be accurate to better 
than 0. 5%. Considerably less accuracy is required as one 
moves away from the maximum mass; it is only 5% for 
the 0.8M Q stars. Given the potential multiple sources of 
numerical error, from the initial data to the numerical 
integrations, it is by no means evident that existing 3 + 1 
numerical relativity codes can set up and track such a 
strong-field collision to 0.5% accuracy JToj] . 

The conjecture put forth in Ref |Q that when two neu- 
tron stars collide head-on from rest at infinity, sufficient 
thermal pressure may be generated to support the hot 
remnant in quasi-static equilibrium against collapse may 
be generally correct, despite the report in Ref of a 
numerical example to the contrary. The simulation cited 
involves the collision of high-mass stars near the maxi- 
mum mass [ pT| , which requires very high numerical ac- 
curacy to simulate reliably In the unreported case 
involving lower-mass stars whose combined mass still ex- 
ceeds the maximum mass, the conjecture appears to be 
corroborated. The accuracy requirements for the lower- 
mass case are far less stringent, and presumably more 
easily attainable. Evaluating the conjecture at high mass 
may have to await more sophisticated code development 
and/or larger machines (l3J|. Even with such advances, 
evaluating the conjecture in this domain may still require 
a careful limiting proceedure based on a sequence of runs 
with increasing initial separation between the initial stars 
to be confident that the initial data correctly corresponds 
to configurations which infall from rest at infinity. 

The 0.8Mq simulation involves stars that have masses 
far below the maximum mass. As a result, this simula- 
tion is likely to be more reliable. If the delayed collapse 
outcome found for this collision is indeed correct, then 
we can already speculate that for neutron stars governed 
by realistic equations of state juj , which currently give 
maximum masses in the range 1.8M Q — 2.3M Q , the colli- 
sion of two realistic (low-mass) 1.4M Q stars from rest at 
infinity likely will lead to delayed collapse in accord with 
the conjecture. 
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1.4M Q . 
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Accounting for these effects places restrictions on the pos- 
sibility of delayed collapse, as discussed in Ref Jl[. 
To test the conjecture, the initial data must not only sat- 
isfy the constraint equations to high precision, as they are 
reported to do in Ref , but must also reliably represent 
the matter and fields at finite separation that arise from 
infall from rest at infinity. The final outcome is increas- 
ingly sensitive to the approximations, like those adopted 
in Ref as the configurations approach the maximum 
mass; see text below. 
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Miller, W. Suen and M. Tobias 1998, |gr-qc/9811015| ). In 
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ficult to measure the total gravitational radiation loss, 
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maximum mass; see || and |n| for discussion. 
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[10] The assessment of Miracle Max, the Miracle Man, upon 
dispensing his latest miracle pill in the movie 
The Princess Bride, comes immediately to mind: 
Valerie [Max's wife]: "Think it will woik [sic: work]?" 
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energy loss is / = AE/M ~ 10~ 3 for the 1.4M Q col- 
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rium solution and delayed collapse is possible ; cf. Eqn. 

<B> 

[12] It is not clear how useful it is to attempt to distinguish, 
as in Ref 0|, the various timescales associated with the 
infall velocity, sound velocity, shock velocity, and ther- 
malization speed in attempting to account for the ob- 
served simulation outcome. To the accuracy that these 
timescales can be defined, they must all be comparable 
on dimensional grounds, given a head-on collision of equi- 
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opment by the authors of Ref p| , may prove necessary. 
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